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ABSTRACT 

We investigate the instantaneous star formation rates (SFR) and extinction properties for a large 
(N = 274), near-infrared (NIR: 2.2/im) + mid-infrared (MIR: 24/xm) selected sample of normal to 
ultra-luminous infrared galaxies (ULIRGs) [10 9 < Lir/Lq < 10 125 ] with (z) ~ 0.8 in the Spitzer 
Extragalactic First Look Survey (FLS). We combine Spitzer MIPS 24/im observations with high- 
resolution, optical Keck Deimos spectroscopy to derive optical emission-line (Ha, H/3, [Oil]) and 
infrared star formation rates (SFR opt & SFRm, respectively). Direct comparison of these SFR 
diagnostics reveals that our sample exhibits a wide range of extinction (1.0 < A v < 4.0 mag). This 
is after removing spectroscopic and IRAC color-selected AGN candidates that account for w 12% of 
the sample. Objects with SFRs of a few solar masses per year have A v values consistent with those 
of normal spirals (A v 1.0 mag). By contrast, LIRGs at z > 1, which make up a large fraction of 
our sample, have SFR w 100 M^yr -1 and a mean A v « 2.5 mag. This translates to a 97% mean 
attenuation of the [Oil] AA3727 forbidden line doublet, with the most extreme sources having as much 
as 99.7% of their [Oil] line flux extinguished by dust. Based on a SFRir/ SFR opt diagnostic, wc 
derive an IR-luminosity-dependent A^ R function [Al R = 0.75 * log(Ljn/ Lq) — 6.35 mag] that we 
use to extinction correct our emission line luminosities. Application of this correction results in a 
correlation between SFRm and SFR opt that has a dispersion of 0.2 dex (Semi-Interquartile Range). 
Investigation of the A v dependence on redshift reveals that for a fixed Lm there is no significant A v 
evolution. Comparisons to previous studies reveal that the mean attenuation of our sample is inter- 
mediate between that of local optical/UV- and radio-selected samples and has a marginally stronger 
Lin dependence. 

Subject headings: galaxies: bulges - galaxies: spirals - galaxies: star bursts - galaxies: absorption 
lines - galaxies: emission lines - galaxies: high-redshift 



1. INTRODUCTION 

Numerous investigations over the past decade have 
been directed at measuring the cosmic star forma- 
tion history (SFH) of the Universe. Observations over 
the full spectrum from radio to X-ray hav e been ex- 
ploited to trace st a r formation rates (se e iKennicuttl 
MOflSt ICoudoulll"992t IChosh Rr. WhT^laioTl for reviews 
of the various diagnostics). The most commonly uti- 
lized have historically been Ha and [Oil] emission 
line and UV continuum flux (eg. iMadau et alJ 119961 
iTresse fc Maddoxl 119981 iHqgg et alJ 119981 lYan et al l 



IGlazebrook et al.ll999UAdelberger fc Steidel200ft 
lErb et al.1 12003(1 . This is largely due to their respec- 
tive accessibility via ground-based observing windows 
for local and distant samples and the fact that they 
are relatively direct tracers of massive star formation. 
Unfortunately, optical and UV diagnostics are highly 
sensitive to dust attenuation. Various approaches have 
been implemented to estimate this reddening. The use 
of Balmer line flux ratios is a direct, but observation- 
ally taxing method that requires high resolution spec- 
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troscopy. In the absence of multiple well measured emis- 
sion lines, color- or magnitude-dependent optical extinc- 
tion c orrections l|Rigopoulou et al.ll2000l iHippelein et alJ 
|2003|) . or the UV slope-extinction relation \0-Ajjy\ de- 
rived for starburst galaxie s (eg. iCalzetti et alJ 119941: 
lAdelberger &: Steidell 120001 have also been adopted. 
Though commonly applied to non-starburst galaxies, the 
latter has been shown t o break down for both more and 
less lumin ous sy stems l)Goldader et al.l 120021 IB ell et all 
2002t lB^l200l . 

By comparison, far-infrared (FIR) and radio star for- 
mation rate (SFR) diagnostics have the advantage that 
they are unaffected by extinction. Their shortcoming, 
however, is that for general populations they have more 
complex relationships to the star formation than opti- 
cal emission line and UV indicators. For instance, the 
typically adopted FIR SFR calibration l|Kennicuttll 1 998(1 
is based on the assumption of infinite optical depth and 
100% reprocessing of massive star UV emission into IR 
flux. This is a reasonable assumption for heavily ex- 
tinguished systems, but breaks down for galaxies with 
moderate attenuation. In spiral galaxies for instance, 
counteracting effects of UV radiation leakage and heating 
from old er evolved stellar populations mu st also be con- 
sidered l|Lonsdale-Persson fc Heloul 119871. Finally, de- 
spite our relatively limited understanding of the decime- 
ter radio continuum, it has served as a powerful proxy 
for the IR flux due to th e tightness of radio- FIR corre- 
lation l|Helou et alJll985l) . Surprisingly, this correlation 
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extends to IR luminosities (« 0.01 L*) at which neither 
the bolometric infrared luminosity (i/ij) no r the radi o 
continuum are reliably tracing star formation ( Bell 2003) . 

Comparative studies of these different diagnostics ex- 
ist for a range of galaxy types and sample selec- 
tions. UV and Ha measurements have been found 
to be ge nerally consistent for local samples of normal 
galaxies llSullivan et alJ l200lt IBell fc Kennicuttl l20?)H 
IBuat et alJ I2002t ISullivan et alJ I2004|) . The scatter 
and, in some cases, the offset between the tracers 
are primarily attributed to a combination of uncer- 
tainties in the extinction correction, s tar/dust geom- 
etry a nd the star formation tim e scale l)Helou fc Bicavl 
ITifflSl iBel 120031 ISullivan et all l200l . The inter- 
comparison of UV/optical to decimete r radio emission 
( Sulli van et alJ200lllAfons o et al.l 2003fl and far-infrared 
IRAS observations llCram et alJl998HDopita et al.l2002t 
iKewlev et al.l 120021 iHopkins et al.| 12003(1 for large local 
samples confirms the importance of accurate UV/optical 
att enuation correctio n s. Recent ISO-b a sed s tudies 
bv iRigopoulou et al] l|2000()lCardiel et all <|2003|) and 
iFlores et all l)2004|) probe out to more distant redshifts 
(z ~ 1) but for admittedly small samples of 12, 7 and 
16 sources, respectively. These pioneering works provide 
the deepest IR-based SFR probes of the distant, dustier 
and more actively starforming universe. Consequently, 
they tend to include more extreme, IR luminous galaxies 
than are present in the local samples. 

Due to the absence of a unified multiwavelength pic- 
ture of the SFH, there is considerable debate about the 
star formation density at high redshift. Most contro- 
versy revolves around issues of sample selection effects, 
SFR calibration uncertainties and dust attenuation cor- 
rections. It has become clear that no single tracer is 
applicable for all galaxies. UV & optical diagnostics ap- 
pear to be well suited for low IR luminosity galaxies that 
do not require significant extinction corrections; whereas 
attenuation-free IR and radio diagnostics provide better 
estimates in dusty, actively starforming systems. In this 
work we study the star formation and extinction prop- 
erties of a large, distant, actively starforming population 
by making a direct comparison of optical emission line 
and IR SFR diagnostics. Specifically, we use mid-IR 
observations from the Spitzer Extragalactic First Look 
Survey and deep Keck optical spectroscopic observations 
to achieve an order of magnitude increase in sample size 
over previous high redshift ISO studies. 

This paper is divided into the following sections. A 
summary of the various observational components and 
their basic analysis is given in §2. Calculations of the 
optical and infrared star formation rates are described 
in §3. Our approach for removing contaminating AGN 
is outlined in §4. The comparison between the optical 
and IR SFRs and our derived extinction corrections are 
discussed in §5. The main points of this work are sum- 
marized in §6. Throughout the paper, we adopt the cos- 
mology of n m = 0.27, fi A = 0.73 and H D = 70 km s" 1 
Mpc -1 . 



The Spitzer Extra-galactic First Look Survey (FLS) 3 
region is a ss 4 sq. deg region centered around 
RA=17H8 m 00 s , Dec=59°30'00". It is chosen to lie 
within the continuous viewing zone (CVZ), have mini- 
mum cirrus and no bright radio sources. Observations of 
this field with each of the 4 IRAC and 3 MIPS imaging 
bands was one of the first science tasks undertaken by 
Spitzer. In addition to IR imaging, numerous ancillary 
datasets including radio, optical and near-IR (NIR) data 
have been taken in this field. Brief descriptions of the 
datasets included in this study are given below. 

2.1. Imaging 

2.1.1. Optical 

Optical i?-band imaging of the FLS was carried out 
using the MOSAIC-I camera on the 4-meter Mayall Tele- 
scope at the Kitt Peak National Observatory on four con- 
secutive nights on UT 2000 May 4-7. A 4 x 2 array of 
SITe 2048 x 4096 CCDs provides a 36' x 36' field of view 
with a 0.258" pixel scale. Tiling 26 individual pointings 
we obtain an i?-band coverage of 9.4 sq. deg, with median 
exposure time per pointing of 1800 sec and typical seeing 
FWHM~ 1.0". The resulting source catalog has a 50% 
completeness limit of R — 24.5 mag (Vega). In addition, 
g 1 & i'-band observations of the FLS were obtained us- 
ing the Large Format Camera (LFC) on the 200-inch 
Palomar Observatory Hale Telescope. These observa- 
tions were taken over multiple observing campaigns from 
August 2001 through June 2004. The total area coverage 
in these bands is roughly 2.0 sq. deg, with comparable 
seeing and resolution to that of the i?-band data. Com- 
pre hensive descr i ptions of these datasets are presented 
in iFadda et all (|2004[) and Glassman et al. (2005, in 
preparation) . 

2.1.2. Near- infrared 

Near-infrared observations were carried out in two 
separate, but complimentary observing campaigns that 
can be characterized as shallow, wide-field and deep, 
narrow-field. In the first, if s -band imaging of a 
1.14 sq. deg region, to a median depth of K s < 19.0 mag 
(Vega) was performed on UT 2001 May 23-26 using the 
Florida Multi-object Imaging Near-IR Grism Observa- 
tional Spectrometer (FLAMINGOS) on the Kitt Peak 
National Observatory 2.1-m telescope. A 2048 x 2048 
HgCdTe Rockwell array provides a 20' x 20' field of view 
with 0.6" pixel scale. Each pointing was comprised of 
50, 30-sec exposures taken in a 25-position dither pat- 
tern. The field was mapped with a 5 x 5 grid pattern 
using half-field offsets (10') between pointings. The me- 
dian exposure time is 2400 sec per pixel, and the median 
stellar PSF over the mosaic is FWHM=1.6". 

In addition to the KPNO dataset, a smaller «45' x 45' 
verification region in the center of the FLS was observed 
with the Wide-field Infrared Camera (WIRC) on the 
Palomar Observatory Hale 200-inch telescope. Observa- 
tions were undertaken over the course of multiple observ- 
ing runs between June 2002 and July 2004. A 2048 x 2048 
Hawaii-II HgCdTe array provides a 8.7' x 8.7' field of view 
with a 0.25" pixel scale. The «0.6 sq. deg area centered 



2. OBSERVATIONS AND REDUCTIONS 



3 For details of the FLS observation plan and the data release, 
see http: / / ssc.spitzer.caltech.edu/fls 
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on the FLS verification region was covered with 34 tiled 
pointings. The average exposure time per pointing is 
3600 sec (120 x 30scc) taken with a 30-position random 
dither pattern, to a median depth of K s < 20.2 mag 
(Vega). A detailed description of all of the NIR obser- 
vations and reductions is presented in Glassman et al. 
(2005, in preparation). 

2.1.3. Mid- Infrared 

The extragalactic component of the Spi tzer FLS is 
comp rised of Infrared Array Camera (IRAC) l|Fazio et al.l 
|2004[) a nd Multiband Im aging Photometer for Spitzer 
(MIPS) (Rick e~et al.l2004^ observations taken in Decem- 
ber 2003 with a total exposure time of 63 hours. The 
MIPS 24^m area coverage was 4.4 sq. deg for the main 
field and 0.26 sq. deg in a deeper verification field, with 
respective 3cr depths of 0.11/0.08 mJy. All data were 
processed and stacked by the data processing pipeline at 
the Spitzer Science Center ( SSC). MIPS photome try was 
performed using Star Finder iDiolaiti et at] <|2000j) . which 
measures profile-fitted fluxes for point sources. A com- 
plete description of the 24^m data reduction and source 
catalog can be found in M arleau et al.l l)2004[) and Fadda 
et al. (2005, in preparation). 

2.2. Spectroscopy 

Optical spectroscopy was obtained with the Deep 
Imaging Multi- Object Spectrograph (DEIMOS; 
iFaber et all 120031) on the W. M. Keck II 10-meter 
telescope. Observations were performed over 3 nights 
from UT 2003 June 27-29. A 1200 line mm" 1 grating 
with central wavelength settings of 7400A and 7699A 
was used with the GG495 blocking filter, resulting in 
a 0.33A pix 1 mean spectral dispersion and a 1.45A 
instrumental resolution. The total spectral range 
observed was 6300 - 9300A; however, the coverage for 
an individual source was limited to 2,630A with a slit 
position dependent starting wavelength. 

A total of 14 multi-slit masks were observed, with 
w 100 1" wide slits per mask. The 5' x 16' slitmasks 
were tiled in 11 unique positions to sample a 25' x 45' 
area centered on the FLS. In Figure ^ mask positions 
are shown on top of the 24yum mosaic for illustration. 
Multiple masks were observed for three of the positions 
in the deepest central region. Table 1 lists the positions, 
PAs and exposure times for the 14 observed masks. 

We require minimum slit lengths of 7" for local sky sub- 
traction and adopt the DEEP2 recommended strategy of 
using tilted slits to better sample and remove sky lines. 
Slit position angles 9 were required to be 10 < \9\ < 25 
degrees from the spatial axis and when possible, were po- 
sitioned along the major axis of elongated galaxies. Of 
the 14 masks, 3 were observed for a total of 10,800 sec 
(3 x 3600sec) and the remaining 10 masks were observed 
for a total of 3,600 sec (3 x 1200sec). Observing condi- 
tions over the course of the run were good, with typical 
seeing FWHM~0.7". For data reduction, we employ the 
DEEP2 spec2d pipeline 4 , which is based on the SDSS 
spectral reduction package. This package performs cos- 
mic ray removal, flat-fielding, co-addition, sky subtrac- 
tion, wavelength calibration and both 2-d and 1-d spec- 
tral extraction. 

4 See http://astron.berkeley.edu/cooper/deep/spec2d/ 




Fig. 1. — MIPS 24/im image of the central 30' X 45' region of the 
First Look Survey Verification region. Positions of 14 5' X 17' Keck 
Deimos optical spectroscopy slitmasks are shown. A summary of 
the spectroscopic observations is given in Table 1. 



TABLE 1 
Spectroscopic Observation Log 



ID 


RA (J2000) 


Dec (J2000) 


P.A. (deg) 


Exposure Time 


Ml 


17:17:08.51 


59:59:53.5 


-39 


3 X 20 min 


M2 


17:16:11.74 


60:00:02.3 


-10 


3 X 20 min 


M3 


17:15:30.86 


60:00:46.0 


-10 


3 X 20 min 


M4 


17:17:07.65 


59:29:58.6 


40 


3 X 20 min 


M5 


17:16:12.00 


59:30:00.0 


40 


3 X 20 min 


M6 


17:15:35.57 


59:29:58.8 


39 


3 X 20 min 


M7 


17:17:41.53 


59:45:26.1 


-10 


5 X 45 min 


M8 


17:17:03.68 


59:45:00.0 


-10 


4 X 45 min 


M9 


17:16:31.14 


59:45:24.8 


-10 


4 X 45 min 


M10 


17:15:53.81 


59:45:05.2 


-10 


3 X 20 min 


Mil 


17:15:12.35 


59:45:00.6 


-10 


3 X 20 min 


M12 


17:17:41.39 


59:45:32.0 


-10 


3 X 20 min 


M13 


17:17:03.68 


59:45:00.0 


-10 


3 X 20 min 


M14 


17:16:31.14 


59:45:24.8 


-10 


3 X 20 min 



2.2.1. Emission-line measurements and Redshift 
Identification 

The 1-d spectral output of the DEEP2 pipeline are an- 
alyzed using an in-house IDL package written by PIC 
& DF. Galaxy redshifts are identified through visual in- 
spection and galaxy template cross correlation. Line flux, 
equivalent width and kinematic measurements are made 
by performing Gaussian profile fits to the 1-d spectra. 
This is done interactively with the user identifying the 
lines to be fit and specifying the spectral region over 
which to perform the line+continuum fit. Our high spec- 
tral resolution (1.45 A instrumental resolution) allows for 
most lines to be modeled individually. One exception is 
the [Oil] doublet AA3726,3729, which is unresolved in 



4 



Choi et al 



~ 1/2 of our sample. Rather than measure each line in- 
dependently, we adopt a double Gaussian profile with a 
fixed line separation of [(1 + z) x 2. 75 A]. We require the 
two lines to have the same FWHM, but allow their flux 
ratio to be a free parameter. 

In the case of Balmer line fits, two-component emis- 
sion-!- absorption profiles are adopted. It is well known 
that nebular Balmer emission lines can suffer from con- 
tamination due to underlying stellar absorption and sub- 
sequently be underestimated. The standard technique to 
correct for stellar absorption is to apply a global Balmer 
line equivalent width correction. Fortunately, our high 
spectral resolution enables us to resolve and directly mea- 
sure the nebular emission and the pressure broadened 
stellar absorption components. This eliminates the need 
for a global correction and results in more accurate line 
fluxes. Figure [21 provides an illustration of two typical 
emission line fits. Figure^ shows a double Gaussian line 
fit to the [Oil] AA3727 doublet, while Figure^ shows an 
emission+absorption two-component fit to H(3. 

2.2.2. Final MIR+NIR Sample 

The spectroscopic observations were designed to target 
a flux-limited NIR-selected sample based primarily on 
our if s -band dataset. Optical g'-, R- & i'-band color in- 
formation was included exclusively to a) clean out stellar 
contaminants and b) apply a rough photometric redshift 
selection, designed to prioritize high redshift sources with 
z > 0.6. The detailed description of this color selection 
will be discussed separately in the spectroscopic catalog 
paper. 

The final sample used for this study includes all galax- 
ies in our spectroscopic sample with a high confidence 
redshift and a significant 24/im detection. In Figure |3| 
the full NIR-selected sample with good spectroscopic red- 
shifts (solid) is shown in contrast to the subsample that 
have 24/im-detected counterparts (dashed). 

The effect of our g' , R,i' , K s color selection is seen in 
both distributions as a sharp dropoff below z < 0.6. Be- 
yond z > 1.1 the impact of our decreasing spectroscopic 
redshift sensitivity is evident. In this regime strong 
sky lines make it challenging to cleanly identify [Oil] 
AA3726,3729. By z > 1.3, this feature has moved out 
of our spectral coverage window for much of our sample. 
The final if s +24/im sample includes 274 galaxies and 
has a median redshift of z = 0.76. 

To investigate how our spectroscopic redshift target se- 
lection is biasing our NIR+MIR distribution, we show in 
Figure 01 the K s versus 24/im flux for the parent photo- 
metric sample (dots) and the targeted spectroscopic sam- 
ple (crosses). Our selection against low redshift sources 
is evident in this distribution by the lack of NIR bright 
(K s < 16.5 mag) spectroscopic targets. 

3. SFR COMPUTATIONS 

3.1. Emission-Line SFR diagnostics 

The Ha A6563 emission line is the most prominent and 
often measured SFR diagnostic in local galaxy surveys. 
Unfortunately, at even moderate redshifts (z > 0.4), this 
line becomes observationally taxing as it moves first into 
the forest of OH sky lines and then beyond the optical 
spectral window. Higher order Balmer lines such as H(3 
and H-f can be related to Ha; however, these are often 



overlooked because their weak line strength and contami- 
nation from underlying stellar absorption make them dif- 
ficult to measure. At redshifts where all the prominent 
Balmer lines become inaccessible (z > 0.9), the [Oil] 
doublet AA3727 can be used based on the Ha/ [Oil] line 
ratio. For this study, we use a combination of all these 
lines to track the SFR over the redshift range of our sam- 
ple. Ha, H(3 and [Oil] arc measured for galaxies in the 
respective redshift ranges z < 0.4, 0.3 > z > 0.9, and 
z > 0.7. When accessible, Hj and H8 are also used as 
secondary diagnostics. In this section, we discuss the flux 
measurement and SFR conversion for each of these lines. 

3.1.1. Line Luminosity Measurement 

The most direct approach for measuring the total emis- 
sion line flux of a galaxy is direct integration of a flux- 
calibrated spectrum. Unfortunately, this is impractical 
for many surveys, given the challenges of properly flux 
calibrating large multi-slit samples. An alternate route 
is to combine line equivalent widths with broadband 

photometry (eg.lHogg. Cohen. Blandford. &: Pahrell998t 

Hopkins. Mil ler. Nichol. Connolly. Bernardi. Gomez. Goto. Tremon 

2003). In our case, we compute fc-corrected absolute 

u',g' & i?-band magnitudes using our observed g' ,R, 

%' and K s photome try and the KCORRECT code 

( Blan ton et alJ 120031) . adapted for our dataset. We 

adopt these restframe magnitudes as approximations 

of the continuum flux density at the rest wavelengths 

of [Oil], HP and Ha (3727A, 4861 A and 6563A, 

respectively). Luminosities of these emission lines are 

then given by: 
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( 3 ) 

where EW is the line equivalent width in angstroms; M 
is the fc-corrected absolute AB magnitude appropriate 
for the rest wavelength of the line being measured; and 
Xrest is the central wavelength of the broadband filter, 
also in angstroms. These derivations assume that the 
continuum flux at a given emission line wavelength is 
well approximated by the flux density at the correspond- 
ing broadband effective wavelength. An additional color 
correction can be made to account for the wavelength 
difference between the line and the filter effective wave- 
length; however, we find this correction to be negligibly 
small (of order a few percent), so it is excluded. 

A major benefit of this approach is that it alleviates 
the need for aperture/slit corrections that can plague di- 
rect flux measurements. These corrections are small for 
compact galaxies, but can be substantial for extended 
sources. Consequently, aperture loss tends to be redshift 
dependent and has the potential to masquerade as evolu- 
tion. By constrast, line equivalent widths and the scaled 
line fluxes described above are fairly insensitive to slit 
loss. This approach does rely on the underlying assump- 
tion that the average line-to-continuum ratio within the 
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Fig. 2. — Line fits of [Oil] and H/3 for two representative objects in our sample. In panel (a), the [Oil] AA3727 doublet is shown. The 
gaussian profiles of the individual lines (dotted) are shown along with their co-addition (dashed). The lines are required to have the same 
FWHM and a fixed 2.75A restframe separation. In panel (b), individual components (dotted) and the coadded H/3 line-fit (dashed) are 
shown. 
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Fig. 3. — Spectroscopic redshift distribution of our full 676 
galaxy sample (solid) and the 274 galaxy, 24/im-detected subsam- 
ple (dashed). 



slit is the same as that outside the slit; however, the same 
assumption is made when correcting for slit loss. 

3.1.2. Balmer Line SFR diagnostics 

The Balmer line star formation rat e diagnosti c s rely on 
the Ha-SFR calibration derived bv iKennicuttJ (Jl998) : 



SFR Ha = 7.90 x \Q- A2 L Ha E Ha (ergs s" 1 ) (4) 

where Lu a is the measured line luminosity and En a is 
the extinction correction factor, measured at the wave- 



length of Ha. This correlation is based on evolutionary 
synthesis models assuming solar abundance, a Salpeter 
IMF with stellar mass limits of 0.1 < M < 1OOM , and 
T e = 10, 000 K case-B recombination. It also assumes 
that the escape fraction of ionizing radiation from the 
observed galaxy is negligible, and therefore that the neb- 
ular emission traces all of the massive star formation. 

The conversion for higher order Balmer lines are de- 
rived based on the intrinsic case-B recombination line 
ratios. As discussed in H2.2.1I these lines tend to be 
less frequently adopted as SFR tracers due the difficul- 
ties of properly c orrecting for stellar absorption; however, 
iKennicuttl l)1992f) has shown that even with the adoption 
of a mean correction, for strong emission line galaxies 
H/3 can be a reliable SFR diagnostic. Incorporating th e 
expected Ha/Hf3 = 2.87 line ratio l|Osterbrocklll989|h 
we obtain the following relation for SFRjjp'- 

SFR H p = 2.75 x 10- 42 L Hp E Hf3 (ergs s' 1 ) (5) 
where the extinction term Eh/3 is measured at H(3. Sim- 
ilar relations for if 7 and HS are derived based on the 
H-f/H/3 = 0.466 and H5/H(3 = 0.256 line ratios. 

3.1.3. [Oil] Forbidden Line SFR diagnostic 

For the most distant z > 0.7 sources in our sample, we 
compute the SFR based o n the [Oil] AA3727 f orbidden 
line doublet, adopting the iKewlev et all l)2004|) calibra- 
tion: 



SFR [OII] = (0.66±0.17)xl0" 



L[oii]E[oii] 



(ergs s 1 ) 

Following the notation above Lrpu] is the measured line 
luminosity and Eton] is the extinction correction fac- 
tor at A3727. An important difference with respect to 



6 



Choi et al 



^ 1000 

x 
3 



3. 

CO 



100 



"i — i — i — r 



"i — i — i — r 



i — i — r 



i — i — i — r 



"i — i — r 



*x x . 



* 



\ ■'*#*■$ *x 



x x * ■■x.tfV >#.xx *■«,.■- 

.■' '.X : / -* S X X ■. 

x x ^ftfe* ** x * 



X 
.X 



<X x x ^ X sf 
< . X x x* ? : .* 

x . x x x x 



■X. x 

X 



J I I L 



J I I L 



J I I L 



J I I L 



J I I L 



15 



16 17 18 

K (Vega) [Mag] 



19 



20 



Fig. 4. — MIPS 24/xm flux versus i\s-band mag for the spectroscopic (crosses) and parent photometric (dots) samples. Comparison of 
the two distributions illustrates that the parent population is being reasonably well represented with our spectroscopic sample for all but 
the brightest sources (K s > 16.5 mag). 



most previous [Oil] calibrations is that Equation does 
not include any assumptions about the differential red- 
dening between Ha and [Oil] of the source . By con- 
trast, the standard conversion of |Kennicu tH II1998D ( cf. 
lOallagher. Hunter, fc Bushousdll989HKennicutti ri992) : 



SFR [OH] = (1.4 ±0.4) x lO- 41 L [OII] E Ha 



(ergs s x ) 
(7) 

requires the calibration of a mean reddening between 
Ha A6563 and [Oil] AA3727. The assumption of an av- 
erage f/a-to-[OII] reddening is reasonable for samples 
with a narrow range of intrinsic extinction; however, it 
would lead to systematic, extinction-dependent errors for 
more general samples with broad E(B — V) distributions. 
Also, since the current sample has a considerably brigh ter 
mean IR luminosity th an that of the|Kcnnicutt (1992) or 
Galla gher et al.l (^89) calibration datasets, using Equa- 
tion [7| would tend to underestimate both the total [Oil] 



extinction and star formation rate. 

It is worth noting that Equation [S] ignores metallic- 
ity effects on the [OII]-to-Balmer line luminosity ratio. 
Kewley+04 does characterize the abundance dependence 
of this ratio; however, due to our lack of metallicity 



measurements, we use their SFR 



[Oil] 



calibration that 



adopts a mean abundance value of the Nearby Field 
Galaxy Survey (NFGS). Though a detailed investigation 
of this issue is beyond the scope of this work, we can 
use their findings to estimate the impact of metal abun- 
dance on our derived SFRs. Over t he metallicity r ange 
8.0 < log{0/H) + 12 < 9.0 Chased o nlMcGaugbll99ll B.™ 
calibration) of the NFGS sample, iKewlev et all l)2004ft 
find that the [Oil] /.Ha ratio exhibits a metallicity de- 
pendent range of ±0.2 dcx. The adoption of a mean 
abundance introduces < 0.08 dex of scatter into their 
[OII]/i/a correlation. We expect a comparable impact 
on the scatter of our SFB[oii] diagnostic. At that level, 
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it does not have a significant effect on our SFR([OII]) 
or extinction uncertainties. 

Finally, in Figure|3]we combine the above emission line 
diagnostics to look at the optically-derived, extinction- 
uncorrected SFR versus redshift. Comparison of the par- 
ent NIR-selected sample with the 24/iin-detected sub- 
sample reveals that latter has a higher mean uncorrected 
SFR opt . This not an unexpected result, since IR flux 
traces star formation. A bit surprising is the degree of 
overlap between the SFR opt distributions of IR-detected 
and non-detected sources. The fact that at any given 
redshift, SFR opt (uncorrected) provides little indication 
of whether a source will be IR luminous or not suggests 
an IR-luminosity dependent attenuation. 

3.2. IR SFR Calculations 

An alternative SFR tracer that is unaffected by ex- 
tinction is the infrared luminosity The IR component of 
galaxy SEDs can be decomposed into three main dust 
emission components: 1) a near blackbody emission pro- 
file of a thermally heated 'cold' big grain (BG) dust com- 
ponent; 2) a near blackbody component of a stochasti- 
cally heated 'warm' very small grain (VSG) dust; and 
3) molecular pol ycyclic aromatic hydrocarbon (PAH) 
emission features ijDesert et al.ll990HDale fc Heloul2f)02l 
hereafter DH02). The primary heat sources for these 
dust components are stellar radiation from young stars, 
an older evolved stellar population and AGN. 

In dusty, high-opacity systems where the dominant 
heat source of the IR dust radiation is young stars (ie. 
Starbursts and LIRGs) the IR luminosity is expected to 
be an excellent tracer of the instantaneous SFR. In these 
situations, the conversion of Ljr to a star formation r ate 
can be made with the calibration of lKennicuttJ 1^998), 

SFR IR = 4.5 x lO- 44 L m (ergs s' 1 ) (8) 

where L\r is defined as the integrated luminosity from 
8 — 1000/im 5 . Equation [H] is based on the assumption of 
solar abundance, a Salpeter IMF (0.1 — 100 M Q ) and an 
optically thick dust distribution. It is consistent (±30%) 
with published calibra tions of comparably selected sam- 
ples ijKenniciIttl 1 1 9981) : however, its extension to more 
general galaxy populations should be made with cau- 
tion. Specifically, the assumption of high optical depth 
places an important limitation on its application to nor- 
mal spiral galaxies. Large UV/optical escape fractions, 
whether due to low dust opacity or dust/star geome- 
try, will cause SFRjr to underestimate the true SFR. 
On the other hand, heating of the diffuse ISM by an 
older background popul ation can contribute significantly 
to the IR l uminosity ijLonsdale-Persson Sz Heloul Il987t 
lHeloull2000|) . resulting in an overestimate of the SFR. At 
the high IR luminosity extreme, the SFRir calibration 
runs into the problem that many ULIRGs derive a signif- 
icant fraction of their bolometric luminosity from AGN. 
In these systems, dust heating from a central AGN can 
be the dominant component of Lir. This contribution 
is difficult to quantify, so it is essential to screen these 
sources. In 2]we discuss our approach for removing AGN 
from our sample. 



3.2.1. Bolometric Correction 

To calculate SFRir, we first compute the bolometric 
IR luminosity from our observed 24/im observations. A 
standard approach for deriving Ljr for lo c al gal axies is 
to use the definition of lSanders fc Mirabell ()1996j) . based 
on IRAS 12, 25, 60 & 100/im luminosities: 

L IR [8 - 1000/im] = 1.8 x 10~ 14 x 10 26 x (9) 
[13.482^(12) + 5.162^(25) + 2.582^(60) + 2^(100)] 

where L ll (X) is in units of LqHz^ 1 and Lir is in units 
of Lq. A comparable Spitzer band transformation ex- 
ists (DH02); however, its application for distant galax- 
ies is hindered by the general dearth of multi-band FIR 
photometry. For instance, the bulk of our FLS sample 
is observed, but undetected with 70 & 160/iin imaging. 
Fortunately, it has been shown from IRAS and ISOCAM 
data that the MIR alone is a reasonab le tracer of Ljr 
ijElbaz et alJ2002HTakeuchi et alJl2005|) . We exploit this 
finding and use our 24/im observations, which corre- 
spond to restframe 10 — 24/im over the redshift range 
of our sample, to derive IR luminosities. 

Rather than implement a single Lj\iiR-to-LiR corre- 
lation, we take the a pproach of using template SEDs 
llCharv FJba7ll200Tl hereafter CE01) to derive L IR [8- 
1000/im] from the 24/im fluxes. CE01 have constructed 
a library of 105 flux-calibrated SEDs that are single- 
valued in Lin and cover the spectral range 0.1 — 1000/im. 
They start with model SEDs matched to a sample of 
local galaxies rang ing from a norma l spiral (M51) to a 
ULIRG (Arp 220) (|Silva et al.lll998fl and combine them 
with MIR ISOCAM spectra and FIR model SEDs. They 
then split and recombine the MIR and FIR components 
of these templates to create a library of SEDs that re- 
produce the observed correlations between the various 
IRAS, ISOCAM and SCUBA IR fluxes of local galaxies. 
We implement the se templates in a manner described by 
Elbaz et al.l l|200 2:) as the 'multi-template' approach. For 
each galaxy in our sample, we shift the template SEDs 
to the redshift of that source. We choose the template 
that most closely reproduces the observed 24/im flux and 
use its 12, 25, 60 and 100/mi integrated fluxes to com- 
pute L/fl[8-1000/im] based on Equation The final 
distribution of Lir for the 2f s +24/im-detected sample 
is shown in Figure El This plot of Lir vs. redshift 
reveals that our sample spans a broad range in both 
Lin and redshift but is strongly concentrated around 
Lir/Lq «2.0 x 10 11 and z «0.8. Also, as a consequence 
of our 24/im flux limit, Lir and redshift are correlated. 
The impact of this correlation on our results will be dis- 
cussed further in t]5.2.3l 

It should noted that the application of this CE01 
SED library relies on the assumption that the luminosity 
trends seen locally are representative of our sample, at 
higher redshift. This has been shown t o be reasonable 
at least out to z ~ 0.8 bv lElbaz et alJ l)2002ft based on 
a comparison of their Radio-MIR vs. MIR-FIR correla- 
tions. In the next section we investigate the uncertainty 
in our derived Lir using an independent family of model 
SEDs. 



5 nb. Different authors have conflicting definitions for Lir and 
Lfir- 
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Fig. 5. — Extinction-uncorrected emission line-derived SFR op t versus redshift for the full 676 galaxy, NIR-selected sample (crosses) and 
the subset of 274 galaxy, i^ s +24^tm-detected sample (circles). The integrated SFR op t distributions of the full (solid) and 24^m-detected 
(dashed) samples show that the mean SFRopt for IR-luminous sources is higher than that of the parent population. The lack of a clean 
separation between the detected and non-detected sources is an expected consequence of a L/^-dependent optical extinction. 



The absolute calibration and the intrinsic uncertainty 
of the mid-IR to Ljr correlation is investigated by com- 
paring our computed IR luminosities to those based on 
model SEDs of DH02. Dale and collaborators generate 
semi-empirical 3-1100/im model SEDs for normal star- 
forming galaxies. These models are created by combin- 
ing three-component (large-grain, very small-grain and 
PAH) dust emission curves based on a power-law dis- 
tribution of the dust mass over heating intensity. In 
contrast to the CE01 family of SEDs, which are based 
on a combination of slightly modified empirical SEDs, 
these are built primarily from theoretical model emis- 
sion curves in which the PAH-dominatcd MIR (3-12/im) 
region is replaced with a modified ISOPHOT spectral 
component. This family of SEDs is single- valued in FIR 
color (/„(60/xm)//i,(100/mi)), indicating that agalaxylR 
SED can be uniquely determined with the measurement 



of this single FIR flux ratio. 

Since we lack FIR colors, rather than try to determine 
the best-fit DH02 template for each galaxy, we compute 
Ltr for the entire family of SEDs, normalized to our 
observed 24^tm fluxes. For this comparison, we adopt 
the conversion: 

Lfg[3 - 1100/oti] = [Ci(z)uL v (2Anm) (10) 
+ C2(z)uL u (70fj,m) + ( 3 {z)vL v (160(j,m)] 

where the coefficients [£i(z), £2 (2)^3(2)] are taken from 
DH02. These SEDs represent the span of star- forming 
galaxy types, so the range of Lf^ provides an estimate 
of the bolometric correction error due to our IR SED 
assumptions. 

In Figured^ top), Lf^ (Lir based on CE01 templates) 
is plotted against the family of Lf^ values (dots). The 



SFR and Optical Extinction in the FLS 



9 



10 



13 



10 



12 







3. 10 11 
o 
o 
o 

1 — I 

I 

00 
OS 

_r 



10 



10 



10' 



"i i i i I i i i r 



~\ i i i xn i i r 



X x xx X 



XX 



X 



X X x ~ x Xx: 



X 



X 

xx 



X 

X x 

X 



X # 



X 



X 



x «f£ x- 
x X x x 



* x 

xx* 



J I I I I I I I L 



X 



J I I L 



J I I L 







0.5 1 
Redshift 



50 



Fig. 6. — Lju versus redshift for our full (N=274) _R" s +24/im-detected sample with accurate spectroscopic redshifts. The Ljr distribution 
is also shown in the right panel. 



bottom panels are the residual plots shown as functions 
of the Ljft {lower-left) and redshift {lower-right). The 
mean ratio of unity for Lf^ j Lf^ indicates that the inde- 
pendently derived CE01 and DH02 Lir values are consis- 
tent on average. The range of this ratio suggests a mean 
LIR uncertainty of «0.2 dex for the sample as a whole, 
and «0.3 dex for the most distant sources at z > 1. At 
the z k0.5 the uncertainty is minimized, indicating that 
the L mir — L[r correlation is tightest for galaxies mea- 
sured at restframe wavelengths of A«15/im. It is worth 
noting that Lf§ [3-1 100/zm] and Lf^[8-1000/xm] are de- 
fined over different wavelength ranges; however, since the 
flux between 3-8/im and 1000-1100/im is of order a few 
percent of the total bolometric luminosity (DH02) this 
difference is ignored here. 

4. AGN CONTAMINATION 

The conversions of both emission line and Lm lumi- 
nosity to a star formation rate hinge on the assumption 



that the dominant ionization and heating source is radia- 
tion from massive young stars. AGN-dominated emission 
line and IR fluxes do not trace star formation and must 
therefore be removed from our sample. Various emission- 
line diagnostics such as [OllV\/H(3 versus [Nil] /.Ha and 
[OIII]/i?/3 versus [Sll]/Ha effectivel y separate popula - 
tions with different ionization sources l|Osterbrocldl 989). 
In cases with more limited spectral coverage, individual 
line ratios suc h as [NelH] / [Oil] have also been success- 
fully utilized l|Kobulnickv fc Kewlevl f2004|) . We do not 
employ these line diagnostics, due to the non-uniform 
rest-frame spectral coverage of our sample. Instead, we 
rely on an IRAC color selection and the visual identifi- 
cation of optical spectral features to identify and remove 
AGN candidates from our sample. 

We first flag sources with obvious AGN signatures such 
as broadened Balmer and/or [Oil] lines or strong Nelll 
or NeV emission. Next, we combine our sam ple with 
IRAC photometry of the FLS l)Lacv et all2005|) to implc- 
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Fig. 7. — Comparison of two independent derivations of the infrared luminosity Lju for sources in our sample. In the top panel, the 
best-fit L^ft for each source in our sample is compared to a range of computed for the full family of DH02 SEDs (black dots). Open 

symbols represent the limiting 'hot' (a =0.06; red triangles) and 'cold' (a =4.00; blue circles) DH02 SEDs. The bottom panels are residual 
plots shown as functions of Ly^ (lower-left) and redshift (lower-right). This figure shows that Ly'^ , adopted as our best-fit bolometric 
luminosity is consistent with the range of Lyy? that one obtains if no Lmir — Lm correlation is assumed. 



ment a 4-band IRAC color selection. It has been shown 
that AGN can be ident ified in the MIR, ba s ed on their 
stron g continuum flux lILaurent et alJl2000t ILacv et aP 
12004: IStern et all 12005(1 . lLacv et all <|2004[) use 4-band 
Spitzer IRAC photometry (3.6, 4.5, 5.8 & 8.0/xm) to iden- 
tify a distinct region in color-color space where quasars 
and AGN are likely to be found. In Figure IBb, we repro- 
duce the IRAC color-color plot from ILacv et alJ l|2004f) 
(Figure 1) for the entire FLS main field. The dashed line 
shows the region expected to be occupied by AGN. In 
Figure IHJ) we show the same plot for our current sample. 
The FLS depth of the IRAC channels 3 & 4 is slightly 
shallower than that of channels 1 & 2, and only 2/3 of 
our sample have clean photometry in all four IRAC bands 
(open circles); for the remaining 1/3 of the sample, limit- 
ing flux ratios are shown (open triangles). Based on the 



IRAC color selection, 9% of each of these subsamples 
fall in the AGN-candidate region. Visual re-inspection 
of the 16 IRAC-selected candidates, reveals that 63% ex- 
hibit some AGN signature in their spectra (6 strong; 4 
moderate). The remaining sources show no obvious in- 
dicator; however, an AGN contribution cannot be ruled 
out, given our limited spectral coverage. We remove all 
16 candidates from the final sample. Investigation of the 
spectroscopically selected candidates shows that 5 of 9 
sources (55%) with strong AGN spectral signatures are 
also IRAC-selected AGN. Ultimately, we find that nei- 
ther the spectroscopic nor IRAC color selection provides 
a complete census of all AGN, so we take the conser- 
vative approach of combining the two methods to clean 
our sample. The impact of this AGN removal is best 
illustrated in the comparison of the two different SFR 
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diagnostics (Fig. |SJ, discussed in the next section. 
5. OPTICAL VS. IR SFR COMPARISON 

Next, we compare the two independently derived SFR 
diagnostics SFR opt and SFRjr described in <J5] for 
the full AGN-cleaned sample. In Figure |§1 we show 
the rcddcning-uncorrected, emission line-derived SFR opt 
versus the IR luminosity-derived SFRir, before and af- 
ter AGN removal. The sample plotted in Figure (left) 
includes the 274 sources with well determined spectro- 
scopic redshifts, emission line and 24/im fluxes. AGN 
candidates are marked based on their spectral (squares) 
or IRAC color identification (circles) (as discussed in 
In Figure M (fight) onr y the final 241 source AGN- 
removed sample is shown. Though the distribution of 
the AGN candidates are not localized to a single region 
of the plot, sources with the most extreme SFR discrep- 
ancies tend to be AGN. This illustrates the potential 
for these sources to be misidentified as heavily obscured 
galaxies. The inclusion of this population would bias the 
sample to appear overly obscured. The fact that AGN 
candidates are found throughout the SFR vs SFR plot 
is not unexpected since the emission line luminosity and 
IR luminosity are loosely correlated in AGN. 

In Figure (right) the uncorrected SFR opt is system- 
atically lower than SFRir . Over the Lir range of our 
sample, the offset spans 2 orders of magnitude and il- 
lustrates the importance of properly accounting for the 
optical extinction. In this section we explore a range 
of different extinction corrections to reconcile the two 
SFRs. In Figure I1UI the ratio of the two diagnostics, 
SFR op t/ SFRjjj, is plotted to better illustrate the sys- 
tematic difference between them. It is evident from Fig- 
ure II Ob . in which no reddening correction has been ap- 
plied that SFR pt underestimates the true SFR by any- 
where from — 2.5 dex. In the remaining panels, we in- 
vestigate the applicability of various fixed extinction cor- 
rections. In Figure HUT b). an extinction of A v — 1.0, rep- 
resentative of normal spiral galaxies is assumed. This is 
a significant improvement over the uncorrected SFR op t , 
especially at low IR luminosity, typical of normal late- 
type galaxies. Beyond Lm > 1O 1O L , however, the two 
SFRs diverge, indicating that this standard correction 
underestimates the attenuation in starburst and brighter 
galaxies. In Figures 110b & d, we adopt more extreme 
corrections of A v = 2.0 mag and A v = 3.0 mag and find 
that neither p rovides a good s olutio n over the full range 
in Lir. The Ca lzetti et alJ 1)2000(1 reddening curve is 
adopted for each of these corrections. We conclude that 
a constant A v correction is insufficient for reconciling 
SFRopt and SFRir. In the following sections, we take 
two independent approaches to characterize the optical 
extinction of galaxies of our sample. In the first, we use 
Balmer decrement and other emission line ratios to com- 
pute extinction, A e ™ . In the second, we adopt SFRjr 
as a proxy for the true total star formation rate and com- 
pute Al R based on the difference between SFR op t and 
SFRir. 

5.1. Optical Extinction Correction I: Balmer Decrement 

A standard method for measuring the internal extinc- 
tion of an individual galaxy is through the comparison 
of observed to predicted Balmer line ratios. We adopt 
this approach for the subset of our sample with multiple 



measured Balmer lines. The color excess E(B — V) gas of 
a source is computed by comparing the observed Balmer 
line ratios {Ff a /F^), (F^/F"?) and (F^/F^) 
with their intrinsic unobscured ratios: 



E(B - V) gas 



2.5 



(F Xl /F X2 ) 

log A/ ,'J (11 



[k(Xi) - k(X 2 )] (Fo 1 / Fo 2 ) 



where (F^/F* 2 ) are the intrinsic unobscured line ra- 
tios based on case -B recombination and T — 10, 000 K 
l)Osterbrocldll989|) . The o bscuration or redden ing curve 
(k(X)) has been derived bv lCalzetti et alJ l)2000|) for star- 
burst galaxies: 



k(\) 



( 2.659(-1.857 + 1.040/ A) + R v 
(0.63^m < A < 2.20/mi) 

2.659(-2.156+1.509/A 
-0.198/A 2 + 0.011/A S ) + R V 
(0.12/xm < A < 0.63/im) 



(12) 



It should be noted that although there are significant 
large-scale differences between the various starburst, 
Milky Way, LMC and SMC reddening curves, over the 
rest-frame spectral region covered by the emission lines 
of in terest (3700A5 , A <6600A), the differences are mi- 
nor (jCalzettil I2001L and references therein). The color 
excess is converted into a wavelength dependent extinc- 
tion based on: 



A(X) = E(B - V) gas k(X) 



(13) 



where A(A) is the mean emission line derived extinc- 
tion in units of magnitudes at the wavelength A. This 
should be distinguished from A 8 (X), the extinction as 
m easured from the stellar continuum. As has been shown 
bv lCalzettH 1)2001(1 these differ by a factor of 0.44 in typ- 
ical local starbursts. 

In addition to using the Balmer decrement, we ex- 
plore the possibility of using the [Oil] /if/3 line ratio to 
compute A e ™ . Though this ratio is less certain than 
that of the Balmer lines, we adopt an empirical value 
for the intrinsic [OH] /if/? line ratio based on the lo- 
cally meas ured reddening-corrected l i ne flu xes. Using 
the NFGS, iKewlev. Geller. fc Jansenl l|2004|) measure a 
mean extinction-corrected [Oil] /H a ratio of 1.2, which 
translates to an [Oil] /H f3 line ratio of: 



F lon ]/Fi 



H0 



= 3.44 



From our full sample, individual reddening corrections, 
A% m , are computed for 24 sources that have multiple 
strong Balmer and/or [Oil] emission lines. Line fluxes 
and SFRopt are then reddening-corrected on a source-by- 
source basis. In Figure ITT1 SFR op t/ SFRjh vs. SFRjr 
for this sample (crosses) is shown before (left) and af- 
ter (right) application of this correction. There is a sig- 
nificant increase in the scatter of the distribution after 
application of the extinction correction. This illustrates 
that due to the enormous uncertainties associated with 
A% m , source- by-source corrections are futile. Despite the 
increased scatter, it is noteworthy that SFR op t becomes 
consistent on average with SFRjr, with the mean offset 
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Fig. 8. — IRAC 4-channel c olor-color p l ots us ed to identify AGN candidates. Panel (a) illustrates how the AGN candidate region of 
color-color space is defined bv lLacv et alj 120041) . All sources with clean 4-color photometry from the f ull 4 sq. deg FLS main field are 
plotted; the region in which quasars and AGN are expected to be found is marked by a dashed line (cf. ISaiina. Lacv. fc Scott 2005, for 
details). Panel (b) shows our sample with the same AGN region marked. The majority of the sample (^ 2/3) has clean photometry in all 
four IRAC bands (open circles). The remaining sources (~l/3) have at least two measured fluxes and therefore at least one limiting color 
(open triangles). Sources that fall in the dashed region of this color-color diagram are flagged as AGN candidates and are removed from 
our final sample. 



reduced from -0.98 dex to -0.07 dex. This indicates that 
in the absence of other extinction indicators, with a large 
enough sample and careful stellar absorption measure- 
ments, Balmer decrement and [Oil]/ Hf3 emission line 
ratios can provide a reasonable first order (Al m ) cor- 
rection for the luminosity range of this sample. The me- 
dian emission line derived extinction of this subsample is 
A e v m = 1.5 ± 1.1 mag with (SFR m ) = 23 ± 8 M Q yr~ x . 
Due to the relatively small size and limited Ljr cover- 
age of this subsample, along with the large uncertain- 
ties, we are not able to derive an L/^-dependent extinc- 
tion correction. Instead, in the next section we adopt 
SFRir as a 'true' SFR and compare it to the extinction- 
uncorrected SFR opt to derive A* R for the full sample. 



5.2. Optical Extinction Correction II: IR vs Optical 

SFRs 

5.2.1. Computing A^ R 

In this section, we take an alternative approach of 
adopting SFR IR as a proxy for the true SFR to deter- 
mine the extinction correction of our sample. Specifi- 
cally, the ratio SFRir over the uncorrected SFR opt is 
used to compute the attenuation of the given emission 
line: 

A(X) = -2.5 * log{SFR opt /SFR IR ) (14) 

where A is the wavelength of the emission line used to 
compute SFR opt . A(X) is converted to a standard visual 
extinction, A v by: 
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Fig. 9. — Comparison of SFR op t and SFRjr before (left) and after (right) AGN removal. In the left panel a sample of 274 sources 
with well-determined emission lines and IR luminosity are plotted as crosses. The 33 AGN candidates identified by their spectral features 
(squares) and IRAC colors (circles) are marked. In the right panel, the final 241 source sample after AGN removal is shown. No extinction 
correction has been applied to SFRopt in either panel. Representative error bars are shown in each panel. The SFR op t error bar is the mean 
la uncertain ty due to the error in line flux measurement. The SFRjr error bar is the mean of the bolometric correction uncertainty, as 
discussed in SI3.2.1I The comparison of these samples illustrates 1) the importance of limiting the AGN contamination and 2) that without 
any reddening correction, SFR op t severely underestimates the true star formation rate, especially for our most IR luminous sources. 



Av ~ k(xj kv (15) 

Extinction derived in this manner may have a large un- 
certainty for individual sources, but this approach should 
produce a reasonable ensemble average. In Figure 1121 
A^ R versus SFRjr and Ljr is shown for our full sam- 
ple. Despite the large scatter in the distribution, a clear 
trend of increasing Al R as a function of SFRir is evi- 
dent. The best-fit line to this distribution is: 

A IR = 0.75(log{S FR m )) + 1.05 (16) 

which in terms of Lir is: 

A IR = 0.75 * log(WL ) - 6.35 (17) 

This fit is limited to starburst and brighter systems 
(Liu > 10 10 L Q ; orange solid line), where SFRjr is 
a good representation of the total SFR. It is extrapo- 
lated to lower luminosities (Lir < 10 10 L Q ; orange dotted 
line). Inclusion of the full sample has only a minor effect, 
slightly steepening the fit. The mean A^ R values (black 
dots) of sources binned in log(SFRm) are overlaid along 
with error bars that show their la dispersion. 



5.2.2. Comparison to Local Samples 

We make a direct compar ison of our A j R fun ction 
to those of local optical /UV iHopk ins et alJ l)2001j) and 
radio-selected samples of lAfonso et all l)2003|) (hereafter 
HOI & A03). We transform their E(B-V) color excess to 
A v based on R v =3.1 and plot them as red dot-dashed & 
blue dashed curves, with dotted segments representing ex- 
trapolations beyond their datasets. For a given SFR, the 
mean A V R of our sample is intermediate between those of 
HOI and A03. The trend in the extinction as a function 
of Ljji is also slightly steeper than the previous rela- 
tionships. To interpret these differences requires a closer 
examination of the various sample selections. 

The HOI SFR dependent extinction relationship is 
based on small (iV»60) local UV+optical selected sam- 
ples (cf. iSullivan et all 1)20011) ). whereas that of A03 is 
derived from a comparably sized radio-selected sample 
that extends to slightly higher rcdshift (z) w0.25 and is 
more actively starforming. The difference in the mean A v 
between the two samples is attributed to the fact that op- 
tical/UV selected samples are biased against heavily ob- 
scured galaxies. Our initial spectroscopic sample being 
NIR-selected (2.2/im), we expect our sample to be signif- 
icantly less affected by obscuration effects than that of 
a UV/optical sample. It is therefore not surprising that 
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Fig. 10. — SFR pt/SFRiR vs. SFRir and L/i{ for the final AGN-cleaned sample, shown for a range of constant extinction corrections. 
In panel (a) no correction for any optical extinction is made. In panel (fc) a uniform standard correction of A v = 1.0 mag, typical of normal 
spiral galaxies, is applied to SFR op t- In panels (c & d) more extreme constant extinction corrections, A v = 2.0 mag and A v = 3.0 mag are 
adopted. Representative error bars are shown in each panel. The SFR op t/SFRjn error bar includes only the mean la uncertainty due to 
the error in line flux measurement, it d oes no t include any uncertainty in SFRir. The SFRm error bar is the mean of the bolometric 
correction uncertainty, as discussed in M3.2.1I With no extinction correction, SFR op t underestimates the true star formation rate by as 
much as 2.5 dex for our most IR luminous sources. Uniform extinction corrections do not adequately represent our sample of normal, 
starburst and LIRG galaxies. 



our mean A(, R is higher than that of HOI. 

The discrepancy between our sample and that of A03 
is more difficult to understand since both radio and 
NIR/MIR selections should be relatively robust against 
obscuration biases. It is suggestive of a selection bias in 
one or both samples. Reliance on optical spectroscopy 
may be one potential culprit. It is possible that the 
dustiest systems are so heavily obscured that either their 
redshifts are indeterminable or their emission lines are 
completely attenuated. We can place some constraints 
on the sizes of these two populations in our sample. Our 
spectroscopic redshift efficiency of 92 ± 5% places an 8% 
limit on the first. This is a conservative estimate since 
some fraction of the 8% are certainly due to our spectro- 



scopic redshift sensitivity function dropping off beyond 
zk1.3. Regarding the second population, 12 sources in 
our sample (< 4%) with identified redshifts, have no mea- 
sured emission line fluxes and therefore only lower limits 
on the derived attenuation. These are shown as arrows 
in Figures l§l fc IT21 and included in our AI r (Lir) esti- 
mate. Deeper optical spectroscopy may reveal these to 
be even more heavily obscured systems. Based on these 
estimates, these two populations are likely to account 
for < 10% of the total population. It is worth noting 
that these biases are not unique to this study since most 
surveys rely on optical spectroscopy for redshift and/or 
emission line diagnostics. 

One potential bias that may overestimate the mean 
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Fig. 11. — SFR pt/SFRjn versus SFRjh and L/^ before and after application of emission line raiio-derived reddening corrections, 
A^ m . A sample of 24 sources is shown (crosses) before (ie/t) and after {right) the application of A% m . For reference, the full sa mple is also 
plotted (small dots). The SFRjr error bars represent the full range of the bolometric correction uncertainty, as discussed in i!3,2,2l The 
SFRopt/SFRjn error bars show the la uncertainty due to the error in line flux measurement, excluding any contribution from SFRir. 
In panel b, individual y-error bars also include the uncertainty in A£ m . Application of AJ m results in a reasonable ensemble- averaged SFR, 
as evidenced by the improvement of (log(SFR op t/ SFRi r)) from -0.98 to -0.07. However, the significant increase in scatter suggests that 
our emission line ratios do not provide optimal corrections for individual sources. 



optical attenuation is due to AGN contamination. Es- 
pecially in the most luminous samples, the contribu- 
tion from AGN becomes increasingly significant. Radio- 
selected, and to a lesser extent, MIR-selected samples 
will tend to have larger fractions of 'active' galaxies than 
comparable depth optical or NIR samples. As seen in 
FigureEl if this population is not sufficiently screened, it 
can result in an overestimate of the sample A v . 

Finally, our A(, R values appear to be consistent with 
tho se derived from an SP SS sample of starforming galax- 
ies l)Hopkins et al.ll2003j) . Though no direct comparison 
is shown here, inspection of their Figure 8 shows that the 
mean A v of their sample is intermediate between that of 
HOI and A03, consistent with the correction derived here 
over a comparable range of SFR. 

5.2.3. Disentangling the Lir vs. Redshift Dependence 

As seen in FigureEl our sample spans a broad range in 
both IR luminosity and redshift; and as a consequence 
of our 24/xm flux limit it exhibits a strong L/jj-redshift 
correlation. So far, we have been assuming that for a 
fixed Lm, Ay is largely redshift independent. Given the 
potential degeneracy between redshift and IR luminos- 
ity dependencies, the validity of this assumption merits 
investigation. 

Though our sample is not ideally suited for a thor- 



ough study of this issue, we can place some limits on the 
redshift dependence by isolating subsets of our data re- 
stricted in IR luminosity and redshift. In Figure ITU we 
take a look at two such slices in Ljr and redshift. The 
first is restricted to 2 x 10 11 < L IR /L Q < 1 x 10 12 , and 
the second to 0.7 < z < 1.0. In the lower panels, plots 
of Ay vs. redshift and Ay vs. Lju for the respective 
subsets reveal that a) for a sample with a narrow range 
of £/_r, there is no strong trend in Ay with redshift; 
and b) over a relatively narrow redshift slice centered at 
z = 0.85, Ay shows a clear correlation with Ljr. Though 
our sample does not allow us to investigate trends over 
the entire redshift or Ljr range of our sample, these two 
slices suggest that to first order the dependence on red- 
shift is small compared to that on Lir. 

6. SUMMARY 

We have combined MIR (24/mi) photometry with high- 
resolution, optical spectroscopy for a large (N = 241), 
Ks + M IPS 24/zm-selected (K s < 20.2 mag; /24 > 
80fiJy) galaxy sample. AGN are removed by implement- 
ing both spectroscopic and MIR color selections. This 
dataset is used to measure the instantaneous star forma- 
tion rate and the mean attenuation of normal through 
IR luminous galaxies (10 9 < Lm < 10 12 - 5 L Q ; (Lm) — 
1.8 x 10 11 J L o ; (SFR) = 31 Afoyr^ 1 ) out to a redshift of 
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Fig. 13. — SFR pt/SFRjn vs. SFRjh and Z//.R of the final AGN-cleaned sample, with various luminosity dependent extinction 
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z < 1.5 ((z) = 0.77). We compare two independent ap- 
proaches of computing the star formation rate. The first 
is based on the IR luminosity, SFRir, and the second is 
based on optical Balmer and [Oil] emission line fluxes, 
SFRopt. 

Comparison of the two SFR diagnostics reveal that 
with no correction for extinction, the optical SFR sys- 
tematically underestimates the IR SFR by as much as 
2.5 dex. This discrepancy has a clear IR luminosity 
dependence that cannot by reconciled with a constant 
A v extinction correction. We take two independent ap- 
proaches to investigate the dust attenuation of our sam- 
ple. 

First, we compute Balmer decrement and emission 
line ratio derived optical extinction A v m on a source- 
by-source basis for a subset of our sample. We find that 



after correction, despite a large scatter in the distribu- 
tion, the ensemble averaged SFR opt is consistent with 
SFRir. The large errors associated with the derived 
A v m , however, illustrate that even with the stellar con- 
tinuum properly measured, the high-order Balmer line 
ratios, such as H(3/Hj and H(3/H5, as well as [Oll]/Hf3 
have only a limited usefulness for extinction corrections 
of individual galaxies at high redshift. This is due to: a) 
the relative weakness of these high order emission lines; 
and b) the limited leverage afforded by the narrow sepa- 
ration of these lines, in comparison to the conventionally 
adopted Ha/H(3 Balmer decrement. 

As an alternative measure of the optical attenuation, 
we use SFR pt/ SFRjr, to derive an IR-luminosity de- 
pendent extinction function, A v = 0.75*log(Lin/ L & ) — 
6.35 mag. In comparing this relationship with local opti- 
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By contrast, the lower-right panel shows that even over a narrow redshift range, Ay has a clear Lm dependence. 



cal/UV and radio-selected samples, we fin d that our Aj R 
functi o n is intermediat e betw een t hose of |Ho^kins et^^l.1 
ll200ll) . ISullivan et all lEfflOl and lAfonso et all il2003ft 
but with a slightly steeper slope. Though it is difficult 
to reconcile these results without a directly overlapping 
sample, we highlight some of th e key differences tha t may 
contribute to the discrepancy. lAfonso et alJ l)2003j) sug- 
gests that the difference between their results is due to 
sele ction bias against hig hly obscured objects in the orig- 
inal Hopkins et al. (2001) sample. Our sample should be 
relatively free o f this selection bias a nd therefore compa- 
rable to that of lAfonso et alJ l|2003fl . The fact that they 
are not is a bit surprising, but may suggest that one or 
both of these samples is still suffering from an additional 
bias. One possibility is that we may be excluding the 
most highly obscured sources due to our spectroscopic 



sensitivities. This is not likely to have a dominant ef- 
fect given the high 92% spectroscopic redshift efficiency 
of our sample and the fact that all spectroscopic samples 
are affected by this bias. Another possibility is that the 
Afons oet alJ l)2003l) sample, due to its radio-selected na- 
ture may contain a larger fraction of galaxies with AGN 
that may be inflating their mean A v . The inclusion of 
even non-dominant AGN may have a similar effect on 
the measured ensemble extinction as that of a heavily 
dust-obscured population. 

We do not attempt to quantify the Ay evolution over 
our full redshift range. However, based on subsamples 
restricted to narrow redshift and Lm ranges, we find no 
evidence for significant Ay evolution. 

The scatter in SFR opt / SFRjh after application of 
our best-fit attenuation correction, is larger than ex- 
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pected from either the line flux errors or the reasonable 
estimates of the IR bolometric correction uncertainty. 
The impact of adopting a mean mctallicity correction 
(< 0.1 dcx) also cannot account for this scatter. This 
indicates that these diagnostics are not providing consis- 
tent tracers of the current SFR, possibly due to differ- 
ences in their dependence on dust geometry and/or SFR 
timescale. 
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